Nonequilibrium Fluctuations, Travelling Waves, and Instabilities in Active Membranes 
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The stability of a flexible fluid membrane containing a dis- 
tribution of mobile, active proteins (e.g., proton pumps) is 
shown to depend on the structure and functional asymmetry 
of the proteins. A stable active membrane is in a nonequilib- 
rium steady state with height fluctuations whose statistical 
properties are governed by the protein activity. Disturbances 
are predicted to travel as waves at sufficiently long wave- 
length, with speed set by the normal velocity of the pumps. 
The unstable case involves a spontaneous, pump-driven un- 
dulation of the membrane, with clumping of the proteins in 
regions of high activity. 

The functioning of active proteins in energy-dissipating 
processes, such as ion transport, protein translocation, 
and biopolymer synthesis, generates forces on the mem- 
branes of the living cell and its organelles ||l|,p| . As the ac- 
tive proteins diffuse around in the membrane, the result- 
ing fluctuations in this force provide a nonthermal source 
of noise for shape fluctuations of the membrane. The 
membranes of a living cell are therefore nonequilibrium 
or active membranes. Although such active, nonequi- 
librium processes are abundant in biological membranes, 
physicists have focussed mainly — with considerable suc- 
cess |^,Q — on the statistical mechanics of membranes at 
thermal equilibrium. There are however reasons MM to 
suspect that nonequilibrium processes are at work even in 
red-blood-cell flicker, traditionally explained as thermal 
equilibrium shape fluctuations 0. The predictions [||J|] 
of fluctuation enhancement in active membranes and the 
micropipette experiments pi3| on membranes laden with 
the photoactive proton pump bacteriorhodopsin (bR) are 
further motivation for our studies. 

In this Letter, we consider the statistical mechanics 
and dynamics of a fluid membrane containing a dis- 
tribution of identical, active pumps free to move in 
the plane of the membrane. By "pump" we mean a 
membrane-spanning protein which, when supplied with 
energy, transfers material (and thus exerts a force on the 
membrane) in one direction only. This ignores the com- 
plexities of some real pumps but is a good description of 
bR. By convention, we shall term the end towards which 
the force acts the head of the pump, and the other end its 
tail. We shall call a protein a "-I- pump" ("— pump") if 
the vector from its tail to its head points parallel to (an- 
tiparallel to) a fixed outward normal of the membrane. 
Flips from -I- to — , prohibitively slow in any real sys- 
tem, are forbidden in our treatment. A membrane will 



FIG. 1. Asymmetric proteins (the black wedges) imbed- 
ded in a fluid membrane are drawn to regions with curvature 
adapted to the protein shape. 



be termed "balanced" ( "unbalanced" ) if the numbers of 
-f- and — pumps are equal (unequal). The unidirectional 
pumping action means that a given pump knows up from 
down. In general, therefore, (i) a pump of a given sign 
will favor one sign of the local mean-curvature H of the 
membrane over the other (Fig. |l|), and (ii) its pumping 
activity will depend on H. Incorporating these effects 
distinguishes the present work from |||,g] . 

Here is a summary of our main results: (i) An ac- 
tive membrane can be either linearly stable or linearly 
unstable against undulation and pump aggregation, de- 
pending on the structural and functional asymmetries 
mentioned in the previous paragraph, (ii) In the sta- 
ble case, for plausible values of the physical parameters 
a "balanced" active membrane has, over an apprecia- 
ble intermediate range of wavenumbers fc, a height vari- 
ance {\hk\'^) oc l/fc**, with a coefficient independent of 
the concentration of pumps. These are truly nonequilib- 
rium fluctuations, with a strength depending on kinetic 
coefficients: (|/ifcP) is proportional to the permeability 
of the membrane. This is broadly consistent with the 
observations of |IG|. (iii) For a stable membrane, the 
longest wavelength disturbances travel as non-dispersive 
waves, i.e., waves with a speed independent of wavevec- 
tor k. The wave speed c is set by the pump activity 
and independent of the membrane elasticity, (iv) In this 
longest wavelength regime, the height variance is much 
smaller, as fc — > 0, than in the equilibrium case; specif- 
ically {\hk\ ) oc p- as k -^ 0. This result depends on 
the effect of nonlinearitics, which also deternrine the fi- 
nal state of the system fl^Jl^ in the unstable case. 

We now construct our model and derive the results 
stated above. We characterize the active membrane by 
its mean curvature H and the areal concentrations v± of 
-f- and — pumps . In a Monge description H ~ — V^/i and 



n±/ yjl + (V/i)^ where h{x,t) is the height of the 
membrane at time t above point x on a two-dimensional 



reference plane, and n±(x, i) are the projections of the 
pump concentrations onto the plane. It is useful to dis- 
tinguish the protein concentration 

n(x, t) = n+(x, t) + n_(x, t) (1) 

and the signed protein density 

777,(x, i) = n_|_(x, i) — n_(x, t) (2) 

with averages mo and no respectively. For convenience, 
we shall work with 



V' 



m 



m 



-,■00 = mo/no. 



(3) 



We begin with a membrane at thermal equilibrium at 
temperature T = /3~^, where the proteins are present 
but not activated. The probability of a configuration 
{/i(x),^(x)} is oc exp(— /3iJ) with a Hamiltonian ||l^ 

E[h,^P] =i I d^x[K{V'^hf + a{Vhf 

+A{iP - tl^of - 2KHi!\7^h] (4) 

to bilinear order in the variables involved. E contains a 
bending energy ||^ with a rigidity k, a surface tension a, 
a compression energy for the signed protein density, with 
osmotic modulus 



Ar^Tno 



(5) 



for no — > 0, and a coupling H which determines the lo- 
cal spontaneous curvature induced |1J] by the presence 
of a protein (Fig. |l]). This last term is one of the con- 
sequences of the directionality of the pumps, and can be 
attributed to a head-tail size difference of magnitude 



no 
We expect £i to be independent of no for no ^ 0. 



(6) 



We now explain why 
tion field (j). Minimizing ( 



1) ignores the protein concentra- 
y) tells us that in thermal equi- 
librium a membrane with a net imbalance tpo of pumps 
will develop a spontaneous curvature 



Hn 



(V^h) 



H 



min E 
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(7) 



Many experiments on the physics of membranes are car- 
ried out on giant vesicles (size ^ 20^m), for which the 
mean curvature and hence, from (Q), ipQ are negligi- 
bly small. Accordingly, we shall work at ■00 = 
for much of this paper. In this limit, the symmetry 
h — > — /i, ijj -^ — i/i (equivalently e —^ —e for each pump 
of sign e = ±) rules out any bilinear coupling -00 in (0) 
for -00 = 0. Thus, in a linearized treatment of a balanced 
membrane, decouples from h and '0. In addition, for 



■00 = 0, the mean normal drift speed of the membrane in 
the active state is zero. Towards the end of this paper 
we shall present some important results for t/jQ ^ 0. 
Equation (0) implies the equilibrium height variance 



{\hk\')^ 



T 



ak^ + Kpffk^ 
at wavevector k, with an effective rigidity 



■^eff 



A 



(8) 



(9) 



independent of sgn(iJ). The dynamics of small fluctua- 
tions is also determined by Heff- We shall assume that 
H is small enough to keep n^ff > [Em, so that the 
presence of the inactive but shape-asymmetric proteins 
merely shifts the value of the bending rigidity. 

In the active state (e.g., when the bR in the experi- 
ments of pn] is illuminated with green light) even static 
quantities such as (|/ifep) must be determined from the 
dynamical properties of the system. To this end, let us 
assemble the ingredients for the equations of motion, to 
leading orders in a gradient expansion, of a membrane 
with active proteins. We continue to work at 0o = 0. 

An isolated active pump with sign e — ± exerts a force 
eFa normal to the membrane, if the local mean curvature 
H = 0. li H ^ 0, the symmetry h — > ~h, e — > — e permits 
an additional contribution £2HFa to the force, where £2, 
whose sign is not fixed by symmetry, is a length char- 
acterizing the sensitivity of the pumping mechanism to 
the bending of the membrane. The force arising from 
the activity of a distribution of + and — pumps with in- 
trinsic concentrations j^±, together with the elastic force 
fei arising from (0), give rise via permeative flow with 
kinetic coefficient /Zp [|8| to a normal velocity 

Vn = tJ-p[{l^+ - V-) + (1^+ + V-)t2H 

+ 0{U+ - U^)'']Fa - fipfel, (10) 

where the 0(i^+ — z^-)"^ (and higher odd order) terms arise 
only if the presence of a given pump affects the activity 
of other pumps. 

Projecting ( |lO| ) normal to the reference horizontal 
plane |1^], defining the natural velocity scale 



Wo = tJ-pFano > 0, 



(11) 



adding the contribution Vhyd due to hydrodynamic flow 
p|,pjj[ and expanding in powers of Vft., we obtain 



dh 
dt 



-^ipdE/Sh + Vhyd + fh 

~ Wo('0 + h'^^h) + Vhyd 

-^ip{-aV'^h + KV*h - kHV^tP) + fh. 



Here the Gaussian zero-mean noise fh has variance pOP] 
2T[iip + {l/4riq)~^] at wavenumber g, and 



(12) 



Vhyd 



d q 
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(13) 



where hq(t) is the spatiaUy Fourier-transformed height 
field and rj the solvent viscosity. The signed protein den- 
sity obeys the conservation law [pl| 






i;^\/h 



1 -I- (V/i)2 



KV^6E/5ip -^ V • fv, 
-V-f0, (14) 



where the vq term comes from projecting (nH) parallel to 
the reference plane, A is the mobility of the pumps, and 
f^ is a spatiotemporally white, isotropic, vector thermal 
noise with variance 2TA. Since the dissipative processes 
that gave rise to fh and f^ are still present, we ignore for 
simplicity possible additional, pumping-dependent bare 
noise sources. In both ( [l2|) and (|lj), the first equality is 
valid for a general energy function E including (pi}} cou- 
plings, while the second (approximate) equality applies 
only to a strictly linearized treatment at -00 = with 
0=1. Equations (|l^ ) and ( |l4|) define our model. 

Before using these equations, it is crucial to note that 
if ^2 or H is sufficiently large and negative, then ( |l^ ) 
and (|lj) predict that a flat membrane with a statisti- 
cally uniform distribution of + and — pumps is linearly 



unstable, even if Kgyj > (see eqn. (g|)), with pertur- 
bations growing at a rate proportional to fc^ at small 
wavenumber k. To see this physically (Fig. 0), imagine a 
region of excess active -I- pumps surrounded by a statisti- 
cally uniform H — mixture, in an initially flat membrane. 
Pumping will pull the -I- region ahead of its surroundings 
leading to mean curvature H = —V^h > around the -t- 
pump. This can lead to an instability in two ways: (i) If 
H < 0, this will attract more + pumps, (ii) if ^2 < 0, the 
activity of the -I- pumps in the region is enhanced. Such 
physical mechanisms could well be involved in processes 
where cell membranes undergo large deformations. Well 
beyond the onset of these instabilities, nonlinearities will 
determine the ultimate fate of the membrane ||ll|,|l2|.^ . 

For the remainder of this paper, we assume parame- 
ter values corresponding to a linearly stable active mem- 
brane. The two eigenmodes for disturbances at small 
wavenumber k are diffusive if fo^i, i'o^2 are small com- 
pared to the pump diffusivity AA, and propagative but 
highly dispersive (wavespeed ~ k) if vq is large enough. 

More important is the the height variance (|ft.fep) at 
wavenumber fc, which is obtained by linearizing and 
Fourier transforming ( [1^ ) and (|lj) in space and time, 
solving for /ik,ij and ipw^uj at wavevector k and frequency 
w, using the statistical properties of the noise sources fh 
and £0 to obtain (|/ik,wP), and integrating over ui. 

For k much smaller than the lesser of 



irjvoH 4r]D^ 

— and kD = . 

A K 



(15) 




FIG. 2. Two mechanisms for instability: (a) If down- 
ward curvature enhances the upward force exerted by a single 
pump, and vice versa; (b) if the curvature produced by pump- 
ing attracts more pumps. 



{\hk\')^ 



T 



FaUolk^ + K^ffk^ 



D^Ik^ 



(16) 



where we have defined ^ = £2 + -^^^i- The thermal wan- 
dering of the membrane is thus suppressed by a pumping- 
induced tension FaUol, but a novel nonequilibrium contri- 
bution to the height fluctuations now appears, mimicking 
a zero-tension equilibrium membrane. Note that the cou- 
pling £ of the pumps to the curvature is crucial here: for 
£ = we would find k~^ behavior as in [g[. For no -^ 0, 
the diffusivity Z?^ = A^ of the pumps approaches that 
of an isolated pump and is hence nonzero and, from (|5|) 
and (|), the length £ ~ ^2 + {n/T)(.i. The coefficient of 
k^ in ( |l6| ) is thus independent of n^ for small tiq. 

If the tension is nonzero, the behavior in ( p^ ) is cut off 
for k less than 



aA 



WIVqUH KKrna 

As a consequence, the excess area |22 






■In- 



(17) 



(18) 



setting the tension a = 0, we find from (12) and (04) that 



can be seen to depend only logarithmically on the pro- 
tein density tiq. Moreover, increasing the solvent viscos- 
ity should decrease fip while leaving D^ unaffected, and 



should thus decrease a. These predictions of (|lg) are 
consistent with the experiments of |lQ|. However, not 
enough is known about the values of parameters such as 
fip, Fa, audi to make a detailed comparison. If instead 
Djp is decreased, {\hk\'^) should increase. Such predic- 
tions of the dependence of static correlations on kinetic 
quantities allow clear tests of our model and of the truly 
nonequilibrium nature of the fluctuations. 

We turn finally to the wavelike modes mentioned at 
the start of the paper. A net excess of pumps of one 
sign (-00 7^ 0) would yield a term vgipo'^'^h in the lin- 
earized (p^). With dl^) this would lead to propagating 
waves with a speed woV'o = Mp^q^^o- Even for -00 = 0, 
the presence of fluctuations means that {4''^) ^ 0. Hence, 
improving upon our strictly linearized treatment by let- 
ting 0^ — > (0^) in ( [lil ) leads to a prediction of waves with 
speed c = voiip^y^"^ = fipFa{m'^y^^, and to a nonequihb- 
rium height variance (|/ifeP) ~ (A*pt^o^)/(^o^2 +-D^)c^fc^. 
It is important to note that these waves only appear for 

wavevectors k < ka^,. = a/ j ^, i for k ^ k^^^^ our ear- 
lier linearized results apply. Hence, for a system with 
small fluctuations (ip^), the linearized results ( [iq ) and 
( |l8| ) could well hold at the experimental wavenumbers k. 

In any case, for 0o 7^ 0; a- complete analysis requires, 
even at a linearized level, the inclusion of fluctuations in 
the protein density n. Replacing "0^ by its average, too, 
is simply a first step in a complete treatment of nonlinear 
fluctuation effects. A detailed consideration of all such 
effects will appear elsewhere |13] . 

We close by summarizing this work. Starting with a 
natural model of a membrane with active proteins, we 
show that such a membrane can be linearly stable or 
unstable to small perturbations in its shape and the dis- 
tribution of proteins. In the stable case, we show that the 
k~^ dependence of the height fluctuations can mimic that 
of an equilibrium membrane, but with an enhanced tem- 
perature, consistent with the observations of pOl. Fur- 
ther experiments, in particular on the dependence of the 
height variance on the protein diffusivity, will provide 
more stringent tests of our predictions. We are currently 
in the process of studying the effects of nonlinearities on 
the scaling of correlation functions in the stable case and 
on the growth in the unstable case ||ll|,|l^,^ . 
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